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1.0  SUMMARY 

This  report  presents  a  numerical  method,  for  solving  the  binary 
diffusion  laminar  boundary-layer  equations.  The  differential  equations  of 
continuity,  momentum,  energy,  and  species  diffusion  are  solved  simul¬ 
taneously  for  two-dimensional  or  axisymmetric  flow.  For  the  case  of 
axisymmetric  flow  the  second-order  transverse  curvature  terms  appearing 
in  the  equations  are  retained. 

Relations  for  mixture  fluid  properties  as  functions  of  local  conditions 
in  the  boundary  layer  are  developed  and  used  in  obtaining  solutions  to  the 
equations.  The  mixtures  considered  in  this  investigation  are  He  -  air, 

Ar  -  air,  and  CO^  -  air. 

Solutions  to  the  equations  are  given  for  flow  of  a  pure  air  free 
stream  over  a  9°  half  angle  cone  at  a  free  stream  Mach  number  of  10 
and  free  stream  Reynolds  number,  based  on  cone  slant  length  of  940.  The 
initial  section  of  the  cone  is  solid  and  mass  injection  is  initiated  at  a 
location  downstream  of  the  tip.  Results  are  presented  and  compared 
for  calculations  both  including  and  neglecting  the  second-order  trans¬ 
verse  curvature  effect. 
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4.  0  NOMENCLATURE 


c. 

l 


f 


g 

H 

h 


K 


1 


K 


2 


k 


L 


Le 


M 

M 


f 

fi 


a 


/  Pe  e 

constant  pressure  specific  heat  of  mixture 
constant  pressure  specific  heat  of  species  f 

constant  volume  specific  heat  of  mixture 
constant  volume  specific  heat  of  species  f 
mass  fraction  of  free  stream  species,  p^/  p 
mass  fraction  of  injected  species,  p^/p 

thermal  coefficient  of  diffusion 

binary  diffusion  coefficient 

stream  function  (equation  16) 
dimensionless  enthalpy  ratio  H/He 
total  specific  enthalpy  of  mixture 
static  specific  enthalpy  of  mixture 
static  specific  enthalpy  of  species  f 

bound  on  trial  solutions  of  momentum  equation 

bound  on  trial  solutions  of  energy  equation 

Boltzmann  constant 

reference  body  length 

Lewis  number  pD,.  C  /X 
fi  p 

molecular  weight  of  species  f 

reduced  molecular  weight  of  species  f  and  i 

count  of  successive  steps  in.  transformed  x  direction 
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P 

Pr 


Rf 


r 

S 

T 

t 

u 

V 


V 


X 


f 


X 


y 

z 

p 

€l 
c  2 
e  3 


fun 

\ 


p 

pf 


pressure 
Prandtl  number 


/,\ 


universal  gas  constant 

distance  from  axis  of  symmetry  to  location  in  the  boundary  layer 

dimensionless  concentration  parameter  (see  equation  27) 
temperature 

transverse  curvature  parameter  (see  equation  25) 
velocity  component  in  x  direction 
velocity  component  in  y  direction 
transformed  velocity  (equation  8) 
mole  fraction  of  species  f 


direction  along  surface  measured  from  leading  edge  or  stagnation 
point 

direction  measured  normal  to  surface 
mass  fraction  ratio  c^/c^ 

e 

dimensionless  velocity  gradient  (see  equation  26) 
convergence  criterion  for  momentum  equation 
convergence  criterion  for  energy  equation 


perturbation  of  initial  condition  on  energy  equation 


transformed  y  coordinate 

maximum  value  of  q  used  in  numerical  solution 
thermal  conductivity  of  mixture 
thermal  conductivity  of  species  f 


viscosity  of  mixture 
viscosity  of  species  f 
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p  mass  density  of  mixture 

4*4  transformed  x  coordinate 

SUBSCRIPTS 

e  evaluated  at  edge  of  boundary  layer 

f  evaluated  for  species  f 

i  evaluated  for  species  i 

n  evaluated  at  station  n 

w  evaluated  at  surface 

m  evaluated  at  free  stream  conditions  or  see  m 

ica 

Primed  quantities  denote  diffe rentation  with  respect  to  r|  or  fj;  barred 
quantities  refer  to  variables  based  upon  transformation  given  by  equations 
(6),  (7),  and  (14),  (15). 
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5.0  INTRODUCTION 

This  report  describes  a  general  method  for  solving  the  laminar 
diffusion  boundary-layer  equations  for  a  binary  gas  mixture.  The  primary 
object  of  this  investigation  has  been  the  development  of  an  exact  method  of 
solution  to  the  governing  flow  equations  applicable  to  arbitrary  two- 
dimensional  or  axisymmetric  shapes  and  a  wide  range  of  flow  conditions, 
such  as  might  be  incurred  during  wind  tunnel  tests  or  actual  flight  condi¬ 
tions. 


Significant  contributions,  which  are  briefly  reviewed  here,  to  the  field 
of  analysis  of  the  binary  diffusion  laminar  boundary-layer  equations  have  been 
made  by  Baron,  Covert,  Culick,  and  co-workers  at  MIT.  Baron  (reference  1) 
developed  the  basic  equations  describing  the  binary  mixture  laminar  boundary 
layer  and  solved  them  for  the  case  of  similar  flows  with  zero  pressure 
gradient,  solutions  were  given  for  Schmidt  number  and  Prandtl  number 
of  unity.  An  integral  method  was  used  to  obtain  the  effect  of  mass  transfer 
downstream  of  the  injection  point.  Cases  were  computed  for  injection  of 
helium  and  carbon  dioxide  into  air.  Covert  (reference  2}  developed  an  approx¬ 
imate  integral  method,  valid  for  small  injection  rates,  which  could  be 
employed  without  the  use  of  large  scale  computers.  Baron  and  Scott  (ref¬ 
erence  3)  extended  the  analysis  to  flows  with  pressure  gradients;  this  study 
was  restricted  to  similar  flows.  Cases  were  computed  for  injection  of 
helium  over  a  temperature  range  for  w'hich  the  vibrational  relaxation  contribu¬ 
tion  to  the  heat  capacity  was  neglected.  A  more  general  analysis  of  the 
problem  (viz.  not  restricted  to  similar  flows)  was  made  by  Moran  and  Scott 
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{reference  4).  The  equations  were  solved  in  terms  of  Crocco  variables 
using  a  finite-difference  technique  developed  by  Fliigge-Lotz  and  Baxter 
{reference  5).  The  method  requires  starting  profiles  and  is  apparently  not 
capable  of  handling  large  changes  taking  place  in  the  streamwise  direction. 
Again  a  helium  -  air  mixture  was  considered  over  a  temperature  range  for 
which  vibrational  relaxation  of  the  air  molecules  could  be  ignored.  Covert's 
method  was  extended  by  Moran  (reference  6)  to  flows  over  a  porous  cone 
with  solid  tip.  In  reference  {7}  Culick  developed  an  integral  method  to 
solve  the  diffusion  equation  for  zero  pressure  gradient  flows.  Velocity 
and  enthalpy  profiles  are  required  for  this  method  which  was  applied  to  flows 
with  helium  injection  into  air.  Subsequently,  Culick  (reference  8)  developed 
a  method  for  zero  pressure  gradient  flows.  This  is  also  an  integral 
method.  It  utilizes  integral  relations  between  heat  and  surface  mass  trans¬ 
fer  and  does  not  require  a  previously  obtained  solution  to  the  energy 
equation.  Results  are  presented  for  both  helium  and  air  injection. 

The  previous  studies  were  primarily  concerned  with  the  application 
of  mass  injection  for  cooling;  moreover,  the  second-order  transverse 
curvature  terms  for  axisymmetric  flow  were  not  included  in  the  equations. 
Recent  experimental  measurements  with  mass  injection  of  helium  and  air 
into  air  (references  9  and  10)  in  hypersonic  flow  over  slender  cones 
indicates  that  mass  injection  has  a  marked  effect  on  drag  coefficients  as 
well  as  heat  transfer,  The  results  reported  in  references  (9)  and  (10)  were 
for  low  Reynolds  number  flows  and  the  total  drag  was  predominantly  skin 
friction  drag.  Hence,  it  is  to  be  expected  that  solutions  to  the  boundary- 
layer  equations  for  flows  with  mass  transfer  will  prove  valuable  for  the 
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prediction  of  the  effects  of  mass  injection  on  drag  coefficients  in  certain 
flight  regimes.  In  addition,  results  reported  in  references  (9),  (10),  and 
(11),  which  were  based  upon,  measurements  for  hypersonic  flow  over 
slender  cones,  indicate  that  the  transverse  curvature  effect  may  be 
significant.  The  theoretical  methods  used  to  treat  this  effect  in  the  above 
references  were  not  exact;  hence,  to  gain  more  insight  into  the  problem, 
the  second-order  transverse  curvature  terms  appearing  in  the  boundary 
layer  equations  have  been  retained  in  the  present  investigation. 

A  numerical  method  based  on  ideas  originated  by  Hartree  and 
Womersley  (reference  12)  is  used  to  solve  the  flow  equations.  The  method 
has  been  developed  for  solution  of  the  incompressible  laminar  boundary- 
layer  equations  (reference  13),  the  compressible  laminar  boundary-layer 
equations  (reference  14),  and  the  laminar  boundary-layer  equations  of  a 
binary  nonequilibrium  dis sociating  gas  (reference  15).  The  primary 
differences  between  the  work  described  in  this  report  and  that  described  in 
reference  (15)  are  the  mathematical  transformations  used  and  the  mixture 
fluid  properties  incorporated  in  the  computer  program.  The  computer 
program  developed  makes  it  possible  to  include  or  neglect  the  second-order 
transverse  curvature  terms  appearing  in  the  equations. 

Numerical  results  for  injection  of  helium,  argon,  and  air  into  air  for 
flow  over  a  9°  half-angle  cone  at  a  Mach  number  of  10  and  free-stream 
Reynolds  number  based  on  slant  length  are  given.  The  free  stream  condi¬ 
tions  are  within  the  range  of  test  conditions  of  the  AEDC-VKF  Low  Density 
Tunnel  L. 
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6.  0  ANALYSIS 
6.  I  Basic  Equations 

The  laminar  boundary-layer  of  a  binary  gas  mixture  with  no  chemical 
reactions  is  considered  for  two-dimensional  or  axisymmetric  flow.  The  actual 
mixtures  studied  are:  He  -  air,  Ar  -  air,  and  CO^  -  air.  Although  thermal 
diffusion  is  neglected  in  the  numerical  calculations  it  is  retained  in  the 
development  of  the  equations.  The  second-order  transverse  curvature  terms 
are  retained  and  handled  through  use  of  a  transformation  introduced  by 
Probstein  and  Elliott  {reference  16.)  The  equations  are  further  transformed 
to  the  form  in  which  they  are  finally  solved  through  use  of  the  Levy  -  Lees 
transformation  (reference  17  ).  The  final  equations  resulting  from  the 
combined  Probstein-Elliott  and  Levy-Lees  transformations  contain  a  common 
term  and  it  is  shown  that  neglect  of  this  term  reduces  the  equations  to  the 
standard  two-dimensional  or  axisymmetric  form. 

The  equations  containing  trans vers e  curvature  in  a  curvilinear  coor¬ 
dinate  system  are  as  follows: 

Overall  Continuity 


h  <rPu>  +  §7  <rPv>  = 0 


(i) 


Momentum 


/  3  u  3  u  \  dp  13/  3u\  .  _ , 

p  (u  ar  +  v  ayr  -di  +  7  a yhd^j 
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Energy 


Species 


»H  ,  8  Hi 

p<u  +  v  bT 


1  8 
r  8  y 


|x  8  H 
Pr  8  y 


+  p  1  -  trju  rc 


+  i  A 

r  8  y 


^CLe-  1)(hf-hi)^+D£T(hf-hi)4f|: 


(3) 


u 


9Cf 

8x 


+  v 


8  cf> 

3y, 


i  JL 

r  8  y 


A 

pDfi  9  y  T  8  yy 


(4) 


6.2  Transformed  Equations 

For  axisymmetric  flow  a  position  in  the  boundary  Layer  is  determined 
by  specifying  the  coordinates  x,  y,  and  r  (see  Figure  1)  where  r  is  given  by 


r  (x,  y)  =  rQ  (x)  +  y  coso- 


(S) 


P  P 


Figure  1.  Coordinates  for  Axially  Symmetric  Body 
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The  body  shape  is  specified  by  the  function  r^  (x).  Equations  (1)  thru 
(4)  can  be  placed  in  an  almost  two-dimensional  form  through  use  of  the 
following  transformation  originated  by  Probstein  and  Elliott  (reference  16). 


ly  =  (r(x,  y)/L)  dy 
ox  =  (rQ  (x)/E)2  dx 


(6) 

(7) 


This  differs  from  the  Mangier  transformation  in  that  r  (x,  y)  appears  in 
equation  (6)  rather  than  r^  (x).  Moreover,  the  Mangier  transformation  is 
normally  applied  only  to  the  boundary  layer  equations  containing  the  first 
order  transverse  curvature  term.  The  y  velocity  component  v  is  eliminated 
through  the  following  relation 


(8) 


Equations  (1}  thru  (4)  in  the  new  coordinate  system  and  in  terms  of 
v  become 


Continuity 


~  (pu)  +  <pv)  =  0  (9) 

8  x  3  y 


Momentum 


P 


(■ 


v 


(10) 
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Energy 


9  H 
9  x 


+  v 


9  H\ 

3y  / 


9H 
9y  + 


H- 


9  u 
9  y 


+  ^  {Le-  1)  (hf 


V 


<hf 


l  3_T 
T  3  y 


(ID 


Species 


/  9cf  Jcf^ 

3 

2 

r 

\(  n 

a  cf  DfT 

9  T  ] 

\U  8*  +  v  Sy) 

=  9  y 

2 

r 

>  Dfi 

3y  T 

3y/J 

o 


The  quantity  (r/rQ)^  is  given  by 


y  cos  a 


(12) 


(13) 


It  can  be  seen  that  neglecting  the  second  term  on  the  right-hand  side 
of  equation  (13)  reduces  equations  (9)  thru  (12)  to  the  two-dimensional  form. 


The  following  transformation  (reference  17)  is  now  introduced  in  order 
to  place  the  equations  in  a  form  more  suitable  for  solution 


%  = 


1 


p  u  u.  dx 
'e  e  n  e 


(14) 


(15) 


and  equations  (9)  and  (  10)  are  combined  through  use  of  a  stream  function 
given  by 

I  „ 

=  ( 2£  )  f  (x.n  )  (16) 
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where 


9  li  —  9  lit 

pu  =  ay  '  pv  55  -  ai 


The  following  relations  are  used  to  express  the  partial  derivative  operators 
in  the  new  coordinate  system 


(Oi  =  (»-) 

vvk 


(b  \  17a  \  /a  s\  /a  \  1 

V>*)y  =  +(a|j- 

V 

(*3\  =  f  (*£\  dy  -  _1  -1°. 

wf/y  Ul>*  J0  \3|/y  2£  {24 


From  equations  (16),  (17),  and  (18)  it  can  be  shown  that 


r,  _/§i\ 

YV: 


f_  u 

x  e 


Introducing  the  stream  function  f,  applying  the  above  operators  and  non- 

_  2 

dimensionalizing  by  the  factor  (4  /pue)  yields 


Momentum 


[C(l+t)  f"]’  + 


f'  )+  ff"  =24 


L  W  _  f  n  U) 
V  34  94/ 
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T| — <x>  : 


(30c) 


(30d) 


(30e) 
(3  Of) 


z 


(30g) 


It  can  be  seen  that  neglect  of  the  quantity  t  in  equations  (22)  -  (24) 
leads  to  the  standard  two-dimensional  form  of  the  laminar  boundary-layer 
equations  in  the  Levy  -  Lees  coordinate  system.  Hence  with  t  set  equal 
to  zero  the  appropriate  transformation  is 


£  = 


I 


2k  . 

p  u  a  r  dx 
e  e  e  o 


Pe  Ue 
(24 


(31) 


(32) 


where  k  =  0  for  two-dimensional  flow  and  k  =  1  for  axisymmetric  flow.  The 
quantity  rQ  does  not  appear  in  equations  (14)  and  (15)  since  it  is  already 
contained  in  the  quantities  x  and  y.  The  second-order  transverse  curvature 
effects  can  now  either  be  accounted  for  or  neglected;  with  t  given  by  equation 
(25)  and  application  of  the  transformation  given  by  equations  (14)  and  (15)  it 
is  accounted  for  and  with  t  equal  to  zero  and  application  of  the  transformations 
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given  by  equations  (3  1}  and  (32)  second-order  transverse  curvature  is  not 
accounted  for.  In  addition,  it  should  be  kept  in  mind  that  when  t  is  zero  and 
equations  (31)  and  (32)  are  applied  the  quantity  r{x,  y)  in  the  original  equa¬ 
tions  (1)  -  (4)  is  replaced  by  r^  (x). 

The  equations  were  programmed  so  that  transverse  curvature  (TVC) 
effects  could  be  included  or  neglected  in  any  given  calculation.  If  TVC  is 
neglected  it  is  understood  that  the  unbarred  variables  apply.  The  expres¬ 
sions  for  certain  quantities  such  as  displacement  thickness  and  the  inverse 
transformation  back  to  the  physical  plane  depend  on  whether  the  barred 
(including  TVC)  or  unbarred  system  (neglecting  TVC)  is  used.  These 
expressions  for  both  systems  are  developed  in  Section  6.5. 

In  order  to  make  the  solution  of  equations  (22)  -  (24)  possible, 
expressions  for  the  fluid  properties  appearing  in  these  equations  are 
necessary. 
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6.  3  Fluid  Properties 

The  fluid  properties  developed  in  this  section  are  those  dealing  with 
a  binary  mixture  of  non-reacting  gases.  The  basic  binary  mixture  is 
obtained  by  the  injection  of  a  foreign  gas,  either  helium,  argon,  or  carbon 
dioxide  into  a  free  stream  consisting  of  air.  The  mixture  properties 
employed  in  this  development  are  functions  of  the  mass  fractions  of  the 
individual  species  and  properties  of  the  pure  species  comprising  the  mix¬ 
ture.  In  a  binary  mixture,  one  of  the  mass  fractions  can  be  eliminated  by 
the  relationship 


c.  +  cr  =  1 

i  f 


(33) 


The  mass  fraction  c{  has  been  eliminated  in  the  present  formulation.  The 
species  equation  variable  z  then  becomes  (as  shown  in  equation  (29)  ) 


z  = 


(34) 


If  the  free  stream  is  uncontaminated,  will  be  one  at  the  edge  of  the 

e 

boundary  layer,  and  if  a  foreign  gas  alone  is  injected  into  the  boundary 
layer,  z  will  be  zero  at  the  wall. 

The  fluid  properties  needed  to  solve  the  laminar  boundary-layer 
equations  of  a  binary  mixture  are:  C  ,  C  ,  q  C  C..  ,  C.r  ,  h, 

P  V  Pi’  Pf'  Vi  ^f 

T 

bf,  hj,  x,  X £>  p.  and  certain  combinations  of  these  that 

give  the  dimensionless  parameters  Pr  and  Le.  Storage  for  the  coefficient 
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of  thermal  diffusion  is  allocated  although  a  numerical  statement  of  it 

is  not  included.  The  fluid  properties  are  valid  in  the  temperature  range 

from  50°K  to  3500°K  (except  for  C  ,  which  is  valid  to  0°K). 

P 

6.  3,  1  Thermodynamic  Properties 


The  individual  species  comprising  the  gas  mixture  are  assumed 
perfect  but  not  necessarily  calorically  perfect.  That  is,  the  specific 
enthalpies  of  the  pure  species  are  not  constant  but  are  taken  as  functions  only 
of  temperature.  Also  the  total  pressure  is  the  sum  of  the  partial  pressures 
of  the  individual  species.  The  expression  for  the  density  p  of  a  mixture  of 
perfect  gases  is  given  by 


P  = 


P 

R77? 


M.M, 

_ i  f _ 

cr(M.-M,)  +  M, 
fir  f 


slugs 


(35) 


Where 


II 

2 

M  .  =  28.  966 

air 

S 

11 

M.  =  39.9d8 

Ar 

=  4.0083 

He 

M  =  44.011 

co2 

The  constant  volume  and  constant  pressure  heat  capacities  of  the 
mixture  are  given  by 


and 


=  (l-cf) 


'V. 


Cf  C 


(36) 


C  =  (  1  -  c,)  C  +  cf  C 

P  f  Pi  f  Pf 


(37) 
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The  constant  volume  heat  capacities  of  the  pure  species  are  obtained  from 
statistical  mechanical  considerations  (reference  18).  The  modes  of  energy 
distribution  considered  are  translation,  rotation,  and  vibration. 


For  the  monatomic  gases,  He  and  Ar,  the  translational  mode  is 
fully  excited  at  room  temperat  ire  to  give  the  value 


C 


V. 

J 


3  RJ_ 
2  M. 


(38) 


and  there  are  no  contributions  due  to  rotation  and  vibration.  For  air  the 
rotational  and  vibrational  modes  must  also  be  used.  The  contribution  of 
electronic  excitations  to  the  heat  capacities  are  negligible  in  the  tempera¬ 
ture  range  of  interest  and  have  been  neglected.  The  constant  pressure 
heat  capacities  of  He,  Ar,  and  air  are  obtained  by  means  of 


C 


(39) 


The  constant  pressure  heat  capacity  for  CO^  was  obtained  from  the  Janaf 
Thermochemical  Data  Tables  (reference  19).  The  C  data  for  air  and 

Pj 

CO^  have  been  fit  to  polynomials  of  the  form 


=  A  +  BT  + 


CT2  + 


DT3  + 


ET  +  FT' 


ft* 


2  oD 
sec  R 


(40) 


The  coefficients  in  equation  (40)  are  given  in  Table  I. 


The  values  of 


for  He  and  Ar  are: 


C 

PHe 


3.  1025  x  10 


ft* 


sec 


R 


(41) 
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3  ft 

C  =  3.  11151  x  10  - - 

Pa  2  °D 

rAr  sec  R 


(42) 


The  static  enthalpy  of  the  mixture  is  given  by 


Where 


and 


(43) 


(44) 


(45) 


Substituting  equations  (40),  (41)  and  (42)  into  (44)  and  (45)  results  in  a 
polynomial  of  the  form 


h. 

J 


AT  +  BT2  +  CT3  +  DT4  +  ET5 


(46) 


for  air  and  CO^  (the  coefficients  for  equation  (46)  are  given  in  Table  II) 
and  for  He  and  Ar  gives 

4  ft2 

hHe  =  3l  1025  x  10  T  JZ  (47) 

sec 


h.  =  3.  1115  x  103T 
Ar  2 

sec 


(48) 


6.  3.  2  Transport  Properties 

The  transport  properties  calculated  in  this  section  are  viscosity, 
thermal  conductivity  and  the  binary  diffusion  coefficient.  The  development 
of  the  thermal  diffusion  coefficient  can  be  found  in  reference  (20). 
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The  viscosity  of  the  mixture  is  obtained  from  Wilke' s  formula 
(reference  21)  applied  to  a  binary  gas  mixture. 


a  = 


^f 


1  +  Gif  3T 

i 


1  +  G 


fi  X, 


(49) 


Where 

(  1  -cf)/M. 

~  c^/Mf  +  (l-c^)/M^ 


(50) 


(51) 


(52) 


and 


are  viscosities  of  the  individual  species 


The  viscosities  of  the  individual  species  are  given  by  the  Chapman-Enskog 
relation  (reference  22). 


j 


0.  625 


kT 


(2,2) 


SI 


(53) 


c? 

For  temperatures  less  than  1000°K,  the  collision  integral  SL'  ’  '  is 
obtained  from  the  Leonard- Jones  (6-12)  potential  as  tabulated  in  reference 
(22);  Sutherland1  s  law,  reference  (23),  is  used  for  air  up  to  a  temperature  of 
1500°K.  For  temperatures  greater  than  1000°K,  the  viscosities  of  He  and 
Ar  are  obtained  from  Amdur  and  Mason  (reference  24);  where  the  collision 
integral  was  evaluated  based  upon  an  exponential  repulsive  potential  of  the 
form 
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_  r/P 

4>  =  Ae  (54) 

The  viscosity  of  CO,,  for  temperatures  greater  than  1000°K  is  taken  from 

reference  (25),  where  a  Morse  potential  of  the  form 

2 

=  A  j^l-e"B(r"c)j  (55) 

is  used  to  evaluate  the  collision  integral.  The  viscosity  of  air  for  tempera¬ 
tures  greater  than  1500°K  is  calculated  by  using  the  mixture  formula 
(equation  (49)  ),  and  taking  the  composition  of  air  to  be  essentially  and 
O^.  The  viscosity  of  is  taken  from  Amdur  and  Mason,  and  the 
viscosity  of  is  taken  from  reference  (26).  The  viscosity  data  for  all  the 
gases  have  been  fit  to  polynomials  of  the  form 

^  =  A  +  BT  +  CT2  +  DT3  +  ET4+FT5^lb  (56) 

The  coefficients  in  equation  (5  6)  are  given  in  Table  III. 

The  thermal  conductivities  of  the  individual  species  are  given  by 
the  Eucken  relation  (reference  27) 


The  above  relation  takes  molecular  vibration  and  rotation  into  account  in 
addition  to  the  translation  contribution  to  the  internal  energy.  The  thermal 
conductivity  of  the  mixture  is  obtained  by  using  equation  (49)  and  replacing 
the  individual  species  viscosities  with  the  conductivities. 
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The  binary  diffusion  coefficients  are  calculated  from  the  first 
Chapman  Enskog  approximation  {reference  22) 

Df.  =  SkV/ltpM^fi*.1*  1J  (58) 


where  is  the  reduced  mass  of  particles  f,  i. 


The  calculation  of  ^  is  made  by  applying  a  reduced  mass 

correction  factor  to  calculated  values  of  D^T  .  (reference  28) 

N2  1 


D  . 
air  i 


(59) 


For  He,  Ar  and  N^,  the  collision  integral  ^  in  equation  (58) 

has  been  evaluated  by  Monchick  (reference  29)  on  the  basis  of  an  exponen¬ 
tial  repulsive  potential  (equation  (54)  ),  where  A  and  p  are  (reference  24). 


A  (ev) 

p(A°) 

He  - 

He 

0.  0386  x  104 

0.  220 

Ar  - 

Ar 

3.  23  x  104 

0.  224 

N2- 

n2 

1.  35  x  104 

0.  263 

To  obtain  the  necessary  combinations  of  these,  a  simple  geometric  mean 
combination  rule  (reference  24)  is  used 
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Using  Monchick'  s  results  and  equation  (58)  values  of  pD^.  were  calculated  and 
fitted  to  a  fifth  degree  polynomial  whose  coefficients  are  given  in  Table  IV. 


For  C0-,  the  collision  integral  was  calculated  based  upon  the  Lennard- 
Jones  (6-12)  potential  (reference  22)  up  to  a  temperature  of  1000°K.  For 
temperatures  greater  than  1000°K,  was  calculated  based  upon  the 
Sutherland  model  (reference  22) 


(62) 


where  D(T)r.  ...  ,  and  Q  are  calculated  from  two  known  tempera- 

fi  -  rigid  sphere 

tures.  pDj.  was  calculated  and  fit  to  equation  (62).  The  coefficients  are 
given  in  Table  IV. 

It  is  stressed  that  the  basic  method  of  solution  is  the  same  for  all 
three  mixtures  considered;  the  only  differences  appearing  in  the  computer 
program  are  the  numerical  constants,  appearing  in  the  previous  equations, 
used  to  describe  properties  of  the  individual  species. 
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TABLE  I 


SPECIFIC  HEAT  AT  CONSTANT  PRESSURE 


?  3  4  5  l 

=  A  +  BT  +  CT  +  DT  +  ET  +  FT  f  — ~ - 

\sec  R 


Gas 

Temperature 
Range  °R 

Air 

0  £  T  <  2000 

2000  <  T  <  6300 

co2 

0  <  T  <  2000 

2000  £  T  <  6300 

■9. 4509125 

,„-4 


■7. 3022675 

,  „-4 


1. 7322782 

,.-6 


■9-  7657438 

,„-10 


1. 7465179 
,„-13 


71993  -4.4825613  7.1435785  -2.4394435  3.4477526  -1.7810174 


17627  7.3287082  -1.342833  1.  3090637  -6.0572879  1.0531063 


1.0499195  -3.4760828  6.1489558  -5.568993  2.033227 
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TABLE  II 


STATIC  ENTHALPY 


h. 

J 


AT  +  BT2  +  CT3  +  DT4  + 


5 

ET  + 


FT 


ro 


Gas 

Temperature 
Range  °R 

A 

R 

C 

D 

E 

F 

Air 

0  £T  <  2000 

6.  0351797 
xlO3 

-4. 7254562 
xl  0"4 

-2. 4340867 
xl  0~4 

4,  3306955 
xlO"7 

-1.  9531487 
xlO"10 

'  ' 

2.  9108631 
xlO'14 

2000  <  T  <  6300 

6.  4371993 
xlO3 

-2. 2412806 

_  1 

xlO  1 

2.  3811928 
xlO-4 

-6. 098608 
xlO'8 

6. 895505 

- 1  2 

xlO 

-2. 968362 

in-16 

xlO 

co2 

0  <  T  ^  2000 

2. 3317627 
xlO2 

3. 6643541 

xlO+  1 

-4.  476109 
xlO"2 

3.  272659 

_  5 
xlO  3 

-1.  2114575 
xlO-8 

1.  755177 

2000  <  T  <  6300 

1. 328997 

4 

xlO 

5. 249597 

-1.  158694 
x!0~  3 

1. 5372389 
xlO"7 

-  1.  1137986 

x  1 0  ” 

3. 388711 

i  n-  16 

xlO 
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TABLE  III 
VISCOSITIES 

(i.  x  10?  =  A  +  BT  +  CT2  +  DT3  +  ET4  +  FT5  (lh  %e - 
J  \  ft  , 

90°R  <  T  <6300  °R 


Gas 

A 

B 

C 

D 

E 

F 

H 

5. 5026909 

7. 9471713 

-2. 2178108 

4.0447782 

-1. 806356 

-1. 01 12767 

xlO'1 

_  3 

xlO 

xlO-6 

,-.-10 

xlO 

xlO~ 14 

xlO'18 

A 

-5. 5369992 

1.0747731 

-4.4382035 

1. 395966 

-2. 0667733 

1.  1450887 

r 

xlO'2 

xlO-2 

xlO"6 

-9 
xlO  7 

in”13 

xlO 

xlO'17 

co.. 

-7. 8193191 

6.7732592 

-17286911 

3.  8700139 

-5. 1304856 

2.  7591624 

L 

xlO'2 

xlO"3 

xlO"6 

in”10 

xlO 

xlO" 14 

.  n- 1 8 
xlO 

Air 

-1.9336805 

8.  9325935 

-3.  918504 

1. 1728501 

-1. 649344 

8.751547 

xlO"2 

xlO'3 

xlO-6 

xlO"9 

in-13 

xlO 

xlO"18 

-  , 
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TABLE  IV 


-9. 7598318 
,.-2 


■3. 2467965 

,.-2 


1. 3094896 

,.-2 


I.  4399797 
,.-3 


4.  8860752 
,.-4 


-5.  6215733 
,.-5 


3.  9911825 


1. 1704382 

.-6 


1. 4178492 

,.-6 


-4. 0516911 


1. 3105574 


3.  8555763 

,.-10 


4. 7053885 
,.-14 


1, 47  35018 
,.-14 


6. 8405177 
,.-14 


-2. 3594221 


-7, 0980313 


-4. 7403394 
,.-18 
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6.4  Overall  Method  of  Solution 

The  method  used  to  solve  equations  (22)  -  (24)  is  based  upon  ideas 
originated  by  Hartree  and  Womersley  (reference  12).  This  method  has  been 
applied  to  the  incompressible  laminar  boundary-layer  equations  (reference 
13),  the  compressible  laminar  boundary-layer  equations  (reference  14),  and 
the  laminar  boundary-layer  equations  of  a  binary  dissociating  gas  (reference 
15).  Many  of  the  details  of  the  method  can  be  found  in  the  aforementioned 
references.  The  essential  differences  between  the  work  reported  here  and 
that  reported  in  reference  (15)  is  in  the  mathematical  transformations  and 
the  fluid  properties  used.  In  addition  a  modification  of  the  integration 
techniques  used  on  the  momentum  and  energy  equations  is  introduced.  The 
derivatives  in  the  transformed  streamwise  direction  appearing  on  the  right- 
hand  side  of  equations  (22)  -  (24)  are  replaced  by  forward  finite  difference 
expressions.  Thus  at  any  given  location  along  the  surface  the  partial 
differential  equations  are  approximated  by  ordinary  differential  equations 
provided  solutions  are  known  at  upstream  locations.  Equations  (22)  -  (24) 
are  rewritten  at  station  n  along  the  body  to  give  the  following  (it  is  to  be 
understood  that  unless  otherwise  specified  quantities  are  evaluated  at  station 
n,  values  at  previous  stations  are  denoted  by  subscripts): 

Momentum 

[C  (1+t)  4>M]  '  =  |3(<K2  +  24>'  +  1  -  (45+r,)  <f>"  +  Rj  (63) 

Energy 

tt1  =  -  (cb  +  T))  41'  +  Rg  (64) 
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Species 


f2 1  =S{<b'+l)  z  -  (<J>+T|)  z'  +  R3 


(65) 


where 


4  '  =  f*  -  1  (66) 

4»  =  g  -  1  (67) 


and 


(68) 


(69) 


The  quantities  4>'  and  4*  are  introduced  to  reduce  round  off  errors  in  the 
numerical  integration  of  equations  (63)  and  (64).  These  values  go  to  zero  at 
the  edge  of  the  boundary  layer  whereas  the  original  dependent  variables 
f’  and  g  go  to  unity.  The  terms  on  the  left-hand  side  of  equations  (63)  -  (65) 
contain  the  molecular  transport  quantities  and  are  composed  of  the  highest 
order  derivatives  in  the  direction  normal  to  the  surface.  The  quantities 
R^,  R^,  and  contain  the  streamwise  derivatives  which  are  first  order 
and  are  expressed  in  finite  difference  form,  for  example 


=  zi 


(4>  *  +  1) 


[a(V 


+  1)  +  b  (4>»  ,  + 

n- 1 


1)  +  c  (4d  + 

n- Z 


ui 


_ 


+  ti)  +  b  (^.j  +  h)  +  c  (4>n_2 


(70) 
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Similar  finite  difference  expressions  are  employed  for  and  R^.  The 
quantities  a,  b,  and  c  depend  on  the  choice  of  step  size  for  the  independent 
variable  £  l  for  the  special  case  of  equally  spaced  steps  of  length  a£ 

a  =  1 .  5  /Af  b  =  -2.0  /A|  c  =  0 .  5  /A  | 

if  solutions  are  known  at  stations  n-1  and  n-2  equations  (63),  (64),  and  (65) 
are  ordinary  in  <j> ,  i|j,  and  z  respectively.  At  station  one,  the  three  point 
finite  difference  expressions  are  replaced  by  two  point  finite  difference 
expressions.  No  difficulty  is  encountered  in  starting  at  station  zero  since 
R^,  R^,  and  R^  are  all  zero  at  £  =0. 

Equations  (6  3)  -  (65)  are  coupled  and  must  be  solved  simultaneously. 

An  iteration  process  is  used  (reference  15).  At  any  given  station  n  the  initial 
fluid  properties  used  are  those  obtained  from  the  final  solution  at  station 
n-1.  The  initial  fluid  properties  are  then  used  to  obtain  a  solution  to  the 
momentum  equation;  the  same  fluid  properties  and  stream  function  obtained 
from  the  solution  of  the  momentum  equation  are  used  to  solve  the  species  and 
energy  equations.  Since  the  fluid  properties  and  stream  function  are 
specified  for  any  given  solution  of  the  species  and  energy  equations,  these 
equations  are  linear  in  z  and  irrespectively.  The  solutions  to  the  species 
and  energy  equations  are  used  to  generate  new  fluid  properties  which  are  in 
turn  used  to  obtain  new  solutions  of  the  species  and  energy  equations.  This 
process  is  continued  until  the  values  of  fluid  properties  converge.  The  same 
stream  function  is  used  for  all  cycles  of  this  iteration  loop.  Upon  convergence 
of  fluid  properties  a  new  stream  function  is  calculated  using  the  latest  fluid 
properties.  This  double  iteration  is  continued  until  convergence  of  solutions 
to  the  momentum  equation  is  obtained.  The  initial  fluid  properties  used  at 
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station  zero  are  obtained  by  assuming  a  linear  variation  between  the  surface 
and  edge  conditions  which  are  known  boundary  conditions. 

6.4.  1  Solution  of  the  Momentum  Equation 

Equation  (63)  is  a  two  point  boundary  value  problem  with  two  conditions 
at  t|  =  0  and  one  condition  at  r|  .  That  is 

d>  <01  =6  (specified) 

w 

4>’  (0)  =  -1 
4>'  <%,)  =  0 

The  boundary  value  problem  is  transformed  to  an  initial  value  problem  by 
hunting  for  the  third  condition  at  r|  =  0,  namely  b”(0)  which  will  give  a 
solution  satisfying  the  condition  <}>'  (q  )  =  0.  With  three  conditions  specified ,4>(0}r 

CD 

4>'(0),  and  <|>  "(0), equation  (63)  can  be  integrated.  The  right-hand  side,  RHS, 
of  equation  (63)  is  first  integrated  and  b"  is  obtained  from 

RHS  dri  +  4>"(0)  (71) 

Successive  integrations  of  b"  then  yield  <$>’  and  d.  It  should  be  kept  in  mind 
that  C  { l+t)  is  known  from  fluid  properties  which  were  obtained  previously. 

Trial  values  of  <£"(0)  are  used  until  three  solutions  have  been  bracketed 
within  the  bound  +  K.  at  p  .  Once  three  solutions  have  been  obtained  and  at 

*  iX> 

least  one  solution  is  either  high  or  low  satisfying  the  condition 

-K.  <b'{r)  )  <K. 

l  »  * 


b 


1 


c  (l+t) 


/ 
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three  point  Lagrangian  interpolation  is  used  to  obtain  the  final  solution 
satisfying  4>'(q  )  =  0.  A  high  solution  is  one  in  which  b'  (q  )>0  and  a  low 

CO  ffi 

solution  is  one  in  which  4>’  (h  )<0.  A  predictor  -  corrector  technique 

CO 

{pp  116-131,  reference  30  }  is  used  to  perform  the  integration  indicated  in 
equation  (71),  Solutions  are  started  at  the  wall  by  means  of  a  Taylor  Series 
expansion  and  a  shorter  step  size  to  maintain  accuracy.  See  references  (14) 
and (20)  for  details  of  the  integration  and  the  Lagrangian  interpolation  processes. 

Under  certain  conditions  it  has  been  found  that  values  of  4>'  (q  )  are 

00 

extremely  sensitive  to  trial  values  of  4)"  (0)  and  difficulty  is  encountered  in 
bracketing  three  solutions  at  q  .  In  previous  work  (references  14  and  15)  the 

03 

original  value  input  for  q  has  been  reduced  when  three  solutions  could  not 

CO 

be  bracketed.  However  in  addition  to  the  explicit  boundary  condition 

p' (q  )  =  0 

CO 

there  exists  the  implicit  condition  that  $  ”  -  0  as  q  »  q  . 

Although  the  solution  is  forced  to  satisfy  the  explicit  condition  the 
implicit  condition  is  satisfied  only  if  q  is  sufficiently  large.  Hence,  q 

ffl  O' 

can  not  be  indiscriminately  reduced  in  order  to  bracket  three  solutions.  A 
method  has  been  devised  to  overcome  this  difficulty.  The  essential  feature 
of  this  method  is  that  once  ^"(O)  bas  been  determined  to  within  a  specified 
accuracy,  the  searching  process  for  no  longer  continues  at  q  =  0  but 
rather  at  some  new  value  of  q.  Solutions  are  generated  by  incrementing 
or  decrementing  <5"  (0),  depending  upon  whether  ^'(q^)  is  low  or  high,  a 
fixed  amount  until  at  least  one  high  and  one  low  solution  is  obtained.  From 
then  on  trial  values  of  4’"  (0)  are  obtained  by  bisection. 
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If  three  solutions  have  not  been  bracketed  within  bound  iK,  and  for 
two  successive  trial  values  (0)  and  <j>."  +  ^(0)  the  following  condition 
exists 


U-'+ 1<°>  -  «»■'«»  I 

0i  +  1<°>  "  4-1  + 


£ 

1 


then  integration  is  begun  at  n  where  q  ,  is  the  maximum  value  of  13  at 

si  s  1 

which  the  last  high  and  low  6’  profiles  differ  by  a  small  amount,  The 

trial  velocity  profiles  A'  agree  at  h  =  0  and  the  difference  between  them 

grows  with  T) ;  hence,  it  can  be  said  that  the  velocity  profile  has  been 

determined  to  within  e  ^  UP  to  qgj.  The  last  three  profiles,  of  which  at 

least  one  is  either  high  or  low,  are  now  used  as  the  basis  for  beginning  the 

search  at  q  Let  the  subscripts  1,  2,  3  refer  to  these  last  three  profiles 

where  1  refers  to  the  low  profile,  2  to  the  middle  profile,  and  3  the  high 

profile.  The  quantities  4*^ (^l s  1 )  and  !,(r|s  ^ )  are  used  as  initial  values  for 

succeeding  integrations  from  q^.  Trial  values  of  4>"{q&^)  are  now  used  to 

obtain  new  solutions.  The  first  new  solution  is  obtained  by  subtracting 

from  or  adding  to  })»  depending  upon  whether  is  high  {>0}  or 

low  (<0),  one  fourth  the  difference  between  4>  ^  (q^  and  <j>'^{qg  ).  The 

process  of  adding  or  subtracting  one  fourth  the  difference  between  the 

last  high  and  low  solutions  to  the  last  trial  value  of  <t>"(qg^)  is  continued 

until  a  high  and  low  solution  are  obtained;  then  bisection  between  the  high 

and  low  is  continued  until  three  values  of  <j>'(q  )  have  been  bracketed  or 

00 

Uj  +  l^sl^  ~  fri’^sl^  1 
<M  +  !<%!>  +  1 
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If  the  latter  condition  is  met  and  three  solutions  are  still  not  bracketed, 

integration  begins  at  q  ?  with  trial  values  of  4>f'(rl  ?)  where  q  ,  is  now  the 

maximum  value  at  which  the  last  high  and  low  solutions  differ  by  This 

procedure  is  continued  J  times;  that  is  until  trial  solutions  have  been 

attempted  from  nsj-  The  quantities  and  J  are  inputs  to  the  program; 

-  8  *6 

typical  values  are  10  ,10  ,  and  4  respectively. 

The  process  of  extending  the  trajectory  of  integration  (ETI)  by 
beginning  from  a  new  location  has  also  been  applied  to  the  species  and 
energy  equations  and  a  schematic  diagram  of  the  process  applicable  to  the 
energy  equation  is  shown  in  Figure  2. 

6.4.2  Solution  of  the  Species  Equation 

With  values  of  the  stream  function  and  fluid  properties  known,  it 
can  be  seen  that  equation  (65)  is  linear  in  z.  The  usual  specified  boundary 
condition  at  the  wall  will  be  z(0).  The  equation  is  second  order  and  a 
value  of  z'(0)  in  addition  to  z(0)  will  provide  the  necessary  information  to 
be  able  to  treat  the  solution  as  an  initial  value  problem.  Two  trial  values 
z’j(0)  and  z2(0)  are  used  to  obtain  and  z^.  These  solutions  are  combined 
linearly  to  give  a  general  solution 

z  =  Az  j  +  Bz2  (72) 

and  then  the  coefficients  A  and  B  are  obtained  from  the  two  boundary 
conditions 

z(0)  =  z  (specified)  and  z(q  )  =  I 

to  give  a  particular  solution. 
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The  quantity  z'{0)  may  be  specified  as  a  boundary  condition  instead 
of  z(0),  In  this  case,  two  trial  values  of  z(0)  are  used  to  obtain  two 
solutions  which  are  combined  to  give  a  general  solution.  The  particular 
solution  is  then  obtained  by  satisfying  the  boundary  conditions 

z'(0)  =  z1  (specified)  and  z{q  )  =  ] 

W  CO 

The  coefficients  A  and  B  in  equation  (7  2)  obtained  from  the  boundary  con¬ 
ditions  are 

1  "  MO 

A  =  - 7 - r  "  > - T  (7  3) 

2,(9  )  -  z  (r,  ) 

J  co  £  CO 

B  =  1  -  A  (74) 

As  with  the  momentum  equation,  the  quantity  R^  is  replaced  by  a 
finite  difference  expression  and  equation  (65)  becomes  ordinary  and  is 
integrated  in  the  rj  direction  at  successive  values  of  The  solution  at 
any  value  of  ^  utilizes  solutions  at  the  two  previous  values  of  £. 

Once  z(0)  and  z'(0)  are  specified  a  solution  z  is  obtained  by  using  the 
specified  initial  conditions  and  integrating  the  right-hand  side  of  equation 
(65)  to  give  U  .  Once  is  determined,  z'  is  obtained  from  equation  (69) 
and  z  is  obtained  by  integration  of  z'.  The  same  integration  formulas 
(pp.  116  -  131,  reference  30)  used  on  the  momentum  equation  are  used  to 
perform  the  necessary  integrations  to  obtain  a  solution  to  the  species  equation. 

6.  4.  3  Solution  to  the  Energy  Equation 

Equation  (64)  is  linear  in  since  fluid  properties  and  the  stream 
function  are  known  from  a  previous  solution.  The  same  general  procedure 
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used  to  solve  the  species  equation  is  used  on  the  energy  equation.  The 
quantity  is  expressed  in  finite  difference  form  to  make  the  equation 

ordinary.  The  specified  surface  boundary  condition  may  be  either  ip(0)  or 
1  ( 0 ) .  If  Lp ( 0 )  Is  specified  two  trial  values  Jjj(O)  and  ^(O)  are  used  to  produce 
solutions  i|/j  and  these  solutions  are  linearly  combined  to  give  a  general 
solution  and  then  the  boundary  conditions 

<M0)  =  4jw  (specified) 

and 

4>(n  }  =  o 

CO 


are  satisfied  to  obtain  the  particular  solution. 

The  general  solution  is  given  by 

4=  A^  +  B^2  (75) 


where  the  coefficients  obtained  from  satisfying  boundary  conditions  are 
given  by 


A 


+l<r,«>  '  +2<n«> 


(76) 


B  =  1  -  A  (77) 

With  given  initial  conditions  ^  (0)  and  d’j(O),  the  right-hand  side  of 
equation  (64)  is  integrated  to  give  it .  Equation  (68)  is  solved  for  ^  j  which 
is  then  integrated  to  give  and  a  similar  procedure  is  used  to  obtain  l[j ^ . 
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If  the  two  trial  solutions  4^  and  are  highly  divergent,  the  linearization 

technique  over  the  complete  range  from  zero  to  q^  can  break  down.  In 

this  case,  a  modified  procedure  is  used  to  obtain  a  solution.  A  set  of 

coefficients  and  Bq  are  computed  as  before  using  equations  (76)  and 

(77);  however,  it  is  assumed  that  the  final  solution  is  valid  only  out  to 

where  q  is  the  minimum  value  at  which  either  ,  or  din  exceeds  the 
PI  PI  1  2 

bound  K2.  (See  Figure  2.) 


That  is 


'l'  =  A  4*1  +  B  4s  for  o 

o  T 1  o  2 

a  new  solution  ^  is  now  generated  from  q 
initial  conditions  where 


£  q  <  q 


PI 


PI 


using  4*(n 


PI 


(78) 


)  and  iKOl  )  as 
t  P1 


iKq  )  = 

PI 


\*.V 


Bo  W 


(79) 


and  )  represents  a  small  perturbation  of  ^'(q  ).  The  new  solution 

PI  PI 

+  3  is  now  linearly  combined  with  ip  to  give  a  continuation  of  the  original 
solution 


=  A1'\*  +  Bj  +  3  for  q  <  q  <  q  (80) 

PI  Pi 

The  process  is  continued  until  a  complete  solution  out  to  ^  has  been 
obtained.  The  quantities  and  are  inputs  to  the  program  and  are 
typically  100  and  10  ^  respectively. 
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6.  5  Boundary  Layer  Parameters 

Once  the  profiles  of  <(>,  z,  ^  and  their  derivatives  have  been  determined 
at  an  x-station,  the  program  determines  the  conventional  boundary-layer 
parameters  of  displacement  thickness,  momentum  thickness,  local  skin 
friction,  and  local  heat  transfer.  The  program  calculates  a  set  of  these 
parameters  for  the  Levy- Lees  transformation  designated  as  t  equal  to  zero, 
and  for  the  Levy- Lees  Probstein-Elliott  transformation  including  second-order 
transverse  curvature,  designated  as  t  not  equal  to  zero  (for  this  case  t  is 
given  by  equation  (25)). 

If  TVC  is  neglected  the  following  expressions  are  used: 

The  displacement  thickness  is  calculated  from  the  formula 


(81) 


when  transformed  to  the  ^  ,  p  -  coordinates  equation  (81)  becomes 


(2£ 


p  u  r 
'e  e  o 


(82) 


The  program  calculates  normalized  displacement  thickness  A*,  which 
is  defined  by 


...  !j  s'" 

A"=(R%)  -f 

r 

o 


(83) 
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where 


Re 

x 


p  U  X 

re  e 


The  program  calculates  y  from  r)  using  the  following  relationship 


Pe  ue 
(ze)1-4 


Yr 

f  -2—  r^  dy 
J  pe  ° 


(84) 


or 


V  =  -till 
y  k 


li 


r  p  u 
o  e  e  ~o 


■n 

/t 


dn 


(85) 


A  normalized  y  is  defined 


Y  =<Rex)^ 


(86) 


and  calculated  from 

Y 


»  f  P*  d 

2k  Ue  )  J  y  dn 

ro  P'e  V  '  o 


(87) 


The  momentum  thickness  is  given  by 


e  = 


7  Kb 


(88) 


or  in  |  >Tj  -  coordinates  by 


8  = 


- k  /f' 

3  u  r  / 
e  e  o  -q 


U-f')  dr| 


(89) 
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A  normalized  momentum  thickness  is  calculated  from 


9 


(Re  ) 
x 


(90) 


The  shear  stress  at  the  wall  is 


r 

w 


(91) 


and  in  £,  coordinates 


r  - 

w 


2  k 

p  u  u  r  f rT 
rw  1  w  e  o  w 


(2£)  ' 


(92) 


Introducing  the  local  skin  friction  coefficient  defined  as 


Cf  = 


w 


J  p  u 
’  r  e  e 


gives 


c 


f 


(93) 


(94) 


The  heat  transfer  at  the  wall  is 


-q 


w 


hi) 


3  c, 


9y 


+  D 


T  1  c?  T 

f  T  By 


w 


(95) 
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-q 


w 


_X  3  h 

C  9  y 
P  7 


(hf 


8Cf  Cr,  T 

1)  ^  _edJ 

3y  \  f 


T 


3  T 


(96) 

w 


In  terms  of  £  ,  ^-coordinates 


r  t  tv 

u.  p  u  H  r 
^  w  w  e  e  o 


w 


Pr  (2£) 

w  3 


lh 


'■■■•  (V) 


w 


Pr  D, 

(Le  -  1)  c  z'  +  — _ — - 
e  ^ 


w 


(97) 


The  Stanton  number  is  defined  as 


St  = 


*w 


p  u  H  (1-g  ) 
e  e  e 


(98) 


(j.  p  r  l 

St=  - — 5 -  )  g  ’ 

Pr  p  (I-g  )  {2£}A|  W 
w  re  w 


(99) 
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If  TVC  is  included  the  following  expressions  are  used: 


The  displacement  thickness  is  given  by 


6*  (l  +  2!"  cos 
o 


¥  cosa¥y 


To  express  6*  in  r\  coordinates 

j—  (ro  y 

=  Vx  +  x 


cose  jdy 


and 


coso- 


(' 


-  1  +  f  1  +  c—a  ( 2  £  ) 

r  p  u 
o  e  e 


>?(Hdi 


A  normalized  displacement  thickness  is  defined  by 


1 

*  =  (Rax>  T- 
o 


and  calculated  from 


&'  = 


(Re  ) 
x 

cos  a 


1  + 


(■•*  I(H  4 


where 


A  =  Zk CQSQ  (2g) 

r  p  u 
o  re  e 


(100} 


(101) 


{102) 


(103) 


(104) 


(105) 
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The  transformation  to  q  and  its  inverse  are: 


P  u 

e 


(106) 


y 


=  {2^) 
p  u 
'e  e 


(107) 


using  equation  (101) 


r 

o 


cose 


(■ 


1+11  + 


ZL,  cosot 

r2  p  u 
o  re  e 


_  54 
(2£  ) 


(108) 


a  normalized  y  is  defined  below 


and  calculated  from 


14 

Y  =  (Re  ) 
x 


Y 


14 

(Re  ) 
x 


COS  Ot 


(109) 


where  A  is  defined  by  equation  (105) 


The  momentum  thickness  is  given  by 


CD 


(HO) 
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and  in  £,  r|  coordinates  becomes 


6  = 


COS& 


-1  + 


(1-  f’)  dq 


fl 


(111) 


A  normalized  momentum  thickness  is  defined  by 


and  obtained  from 


e 


(Re  ) 
x' 


(112) 


(Re  )' 
x 

coscv 


7 

- 

f  r 

(!-£■)  dif* 

(113) 


The  shear  stress  obtained  from  equation  (91) 
in  q  coordinates  is  given  by 


w 


p  fi  U2  r  f" 
rw  w  e  o  w 


-  i 


L  (2£) 


i 


(114) 


Introducing  c 


f 


defined  by  equation  (93)  gives 


c 


f  = 


f" 

w 


(115) 


The  surface  heat  transfer  obtained  from  equation  (95)  in 
t,  9  coordinates  is  given  by 
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■qw  = 


u,  p  u  H  r 
^  w  rw  e  e  o 

LPr  (2|)^ 
w  ^ 


g'  + 


fh,  -  h. 

f  1 


H 


e  /  w 


'Le-1)  c 


f 


T 

Pr  Df 

z»  +  ~  T- 


w 


016) 


Stanton  number  is  given  by  equation  (98)  and  becomes 


St  = 


(i.  p  r 
^  w  w  o 


P'.„  LP.  ('-Bw*  <2I>* 


w 


w 


(117) 
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7.0  RESULTS  AND  DISCUSSION 

The  laminar  boundary-layer  flow  field  was  calculated  for  a  9°  half¬ 
angle  cone  at  a  free  stream  Mach  number  of  10  and  free  stream  Reynolds 
number,  based  on  cone  slant  length,  of  940.  Cases  were  computed  with  and 
without  mass  injection.  The  following  edge  conditions  were  used  for  all 
cases  computed:* 


u  =  7897.0  ft/sec  p  =  3. 485  x  10“ 7  lb  sec2/ft4 

e  e 

H  =  3.  395  x  107  ft2/sec2  T  =  460.09  °R 

e  e 

P  =  0.  2749  lb/ft2  P  =  3.  368  x  10~7  lb  sec/ft2 

e  e 

cf  =1.0 
e 

The  conditions  at  the  wall  are  specified  bv  f  ,  g  ,  and  z  .  At 

w  °w  w 

locations  where  the  wall  is  solid  f  is  zero;  at  locations  where  mass 

w 

injection  takes  place  f  is  given  by  equation  (30a),  If  TVC  (transverse 

w 

curvature)  is  neglected  f  _  is  given  by  a  slightly  modified  form  of  equation 

^■r 

(30a),  shown  below 


f 

w 


1 

<z£V’ 


i 


p  v  r  cbc 
w  w  o 


*  These  values  are  within  the  range  of  test  conditions  of  the  AEDC-VKF 
Low  Density  Tunnel  L. 


i 
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It  should  be  kept  in  mind  that  without  TVC  the  Probstein-Elliott  transforma¬ 
tion  is  not  used  and  the  independent  variables  are  unbarred.  The  value  of 
gw  used  is  based  on  a  surface  temperature  of  1116.  6°R;  for  pure  air  this 
corresponds  to  a  value  for  of  0.2,  The  value  of  is  unity  at  surface 
locations  with  no  mass  transfer  and  is  zero  at  locations  where  there  is 
injection  of  a  foreign  gas. 

The  body  used  for  computational  purposes  is  based  on  a  test  model 
employed  at  AEDC-VKF.  It  is  a  9°  half-angle  cone  with  solid  front  section 
followed  by  a  porous  section  used  for  mass  injection.  Taking  the  slant 
length  of  the  cone  as  l  the  porous  section  begins  at  an  -x.il  of  0.4218.  The 
actual  length  of  the  test  model  is  2.  39?  in.  and  the  free  stream  Reynolds 
number  used  is  based  on  this  length.  In  all  cases  calculated  with  blowing 
the  mass  injection  rate  is  assumed  uniform  from  0.4218/  to  the  base  of  the 
cone  (see  Figure  3). 


Figure  3.  -  Diagram  of  Cone  Model 
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The  results  presented  in  this  section  are  for  calculations  both  with  and 
without  the  TVC  term  included  in  the  equations.  Both  sets  of  results  are 
exhibited  in  order  to  provide  a  basis  for  comparison  of  predicted  values 
obtained  from  a  standard  Mangier  treatment  with  values  obtained  from  an  exact 
treatment  of  TVC. 

The  TVC  term  t  as  can  be  seen  from  equation  {25)  is  of  the  form 


where  the  function  F(|)  is  well  behaved  if  the  radius  of  the  body  is  finite. 
However,  if  the  radius  is  zero,  which  is  the  case  at  the  tip  of  a  cone,  F(|) 
becomes  singular.  Trial  solutions  indicated  that  downstream  solutions  are 
insensitive  to  the  value  of  F(0}.  Numerical  solutions  were  obtained  first  by 
neglecting  F(0)  and  using  various  distributions  of  £  near  the  tip  of  the  cone. 
The  different  solutions  converged  to  the  same  values  downstream  of  the  tip 
and  it  was  shown  that  the  region  of  convergence  moved  upstream  as  locations 
closer  to  the  tip  were  used.  Neglecting  F(0),  however,  did  cause  the  results 
in  a  region  near  the  tip  to  be  of  questionable  accuracy;  for  example,  the  skin 
friction  coefficient  rather  than  decreasing  monotonically  exhibited  a  slight 
dip  and  rise  near  the  tip  prior  to  decreasing.  This  phenomenon  was 
eliminated  by  employing  a  value  of  F(0)  which  is  higher  than  F(^)  where 
is  the  first  location  past  the  tip  at  which  a  solution  is  obtained.  In  all  calcu¬ 
lations  involving  TVC,  the  first  station  past  the  tip  was  taken  at  a  value  of 
£  of  0.  236231  x  10“^;  this  corresponds  to  an  x/£  of  2.  5  x  10"^.  As  the 
solution  proceeded  downstream,  larger  increments  of  £  could  be  handled 
and  a  total  of  27  stations  were  used  for  the  complete  cone.  The  value  of 
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F(0.  236231  x  10"^)  for  the  particular  flow  studied  here  is  16.  1496  and  F(0) 
was  arbitrarily  chosen  as  32.  An  examination  of  equations  (22),  (23)  and 
(24)  indicates  that  as  larger  values  of  F{0)  are  used  the  values  of  f"(0),  g1  (0), 
and  z’  (0)  increase  and  in  the  limit  as  F(0)  »  co  f',  8>  an<^  2  are  unity  for  all 

q  greater  than  zero.  The  use  of  a  large  value  for  F(0)  will  provide  steep 
gradients  at  the  wall  and  hence  provide  reasonable  starting  profiles. 

The  Hartree  Womersley  method  was  originally  developed  to  handle  a 
one  component  gas  (references  14  and  31).  A  different  transformation  than 
is  employed  in  the  present  investigation  was  used  and  the  method  was  pro¬ 
grammed  for  pure  air.  In  the  present  investigation,  if  z  is  unity  and  z'  is 
zero  for  the  whole  boundary  layer,  then  the  binary  mixture  equations  reduce 
to  those  describing  pure  air.  Although  the  solution  of  the  species  equation 
would  yield  a  z  of  unity  throughout  the  boundary  layer  for  a  single  component, 
an  option  was  incorporated  in  the  computer  program  which  can  be  used  to 
bypass  the  solution  of  the  species  equation  and  automatically  yield  a  value  of 
unity  for  z  and  zero  for  z1.  The  fluid  properties  computed  are  then  those  of 
pure  air.  Using  this  option,  the  boundary  layer  with  mass  injection  of  air  on 
the  body,  shown  in  Figure  3,  was  computed  neglecting  TVC  and  compared  to 
results  using  the  method  described  in  references  14  and  31.  The  fluid  property 
relations  used  in  the  two  methods  and  the  mathematical  transformations 
differ;  values  of  displacement  thickness  and  Stanton  number  divided  by 
skin  friction  coefficient  shown  in  Table  V  are  in  fair  agreement.  This 
comparison  was  made  primarily  for  the  purpose  of  checking  the  system 
of  reducing  the  binary  mixture  equations  to  those  for  a  single  component. 
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Figures  4  through  11  and  Tables  VI  and  VII  exhibit  results  from  calcu¬ 
lations  both  neglecting  and  including  TVC  and  employing  helium,  air,  argon 
and  carbon  dioxide  as  the  injected  gas. 

Since  viscous  interaction  is  neglected  in  the  present  calculations,  the 
quantity  (3  appearing  in  equation  (22)  is  zero  and  since  the  free  stream  is 
pure  air  for  all  locations  along  the  cone,  the  quantity  S  appearing  in  equation 
(24)  is  zero.  The  circled  points  in  Figures  5(b)  and  10(b)  indicate  locations 
on  the  cone  at  which  solutions  to  equations  (22)  through  (24)  were  obtained. 
Intervals  of  0,  05  of  the  transformed  y  coordinate  were  used  with  smaller 
step  sizes  (see  references  20  and  31)  near  the  wall  to  maintain  accuracy. 

Probstein  and  Elliott ( refer ence  16}  point  out  that  the  additional  term  in 
the  differential  equations  arising  from  the  inclusion  of  the  transverse  curva¬ 
ture  effect  causes  behavior  similar  to  that  of  a  favorable  pressure  gradient. 
Examination  of  the  figures  and  tables  indicate  that  inclusion  of  the  TVC  term 
causes  predictions  of  skin  friction  coefficient  and  Stanton  number  to  be  higher 
and  predictions  of  displacement  thickness  lower  than  predictions  based  upon 
calculations  which  neglect  TVC;  hence,  the  present  results  agree  with  the 
aforementioned  observation  of  Probstein  and  Elliott. 

The  computed  results  show  that,  in  general,  the  effect  of  TVC  increases 
with  decreasing  molecular  weight;  hence,  its  effect  is  greater  for  helium 
injection  than  injection  of  any  of  the  other  gases.  This  result  is  explained 
by  the  fact  that  the  TVC  term  t  appearing  in  the  equations  contains  the 
quantity  ^p^/pdq  which  increases  with  decreasing  molecular  weight.  More¬ 
over,  the  results  indicate  that  inclusion  of  TVC  predicts  values  of  the  various 
boundary  layer  parameters  which  differ  significantly  from  values  obtained 
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from  calculations  which  neglect  the  TVC  effect.  For  example  (see  Figures 
5(a)  and  9(a)  )  with  no  mass  injection,  the  skin  friction  coefficient  at 
x/i  =  1.0  is  about  50%  greater  when  TVC  is  included.  It  is  also  shown  (see 
Tables  VI  and  VII),  that  the  relative  effect  of  helium  injection  in  comparison 
to  air  injection  is  far  less  pronounced  when  TVC  is  taken  into  account. 

Of  interest  in  connection  with  the  results  obtained  here  are  experimental 
results  reported  by  King  and  Talbot  (reference  9)  which  show  that  for  a  given 
injection  rate  air  has  a  greater  effect  on  the  drag  coefficient  than  does  helium. 
The  present  results  in  this  respect  do  not  agree  with  experiment;  however, 
it  can  be  said  that  the  inclusion  of  TVC  in  the  calculations  greatly  improves 
matters  regarding  agreement  between  theory  and  experiment. 

It  is  felt  that  inclusion  of  the  boundary  layer  displacement  effect  on  the 
external  flow  would  be  one  further  refinement  which  might  clear  up  the 
discrepancy  between  theory  and  experiment.  This  effect  would  introduce  a 
favorable  gradient  and  as  with  TVC  the  effect  would  increase  with  decreasing 
molecular  weight  since  the  pressure  gradient  term  appearing  in  the  momentum 
equation  (i.  e.  ,  see  equation  (22)  )  contains  the  quantity  pe/p. 

The  present  results  indicate  the  importance  of  solutions  utilizing 
accurate  fluid  properties  to  ascertain  the  relative  effects  of  various  injection 
gases  on  the  boundary  layer.  For  example,  a  simplified  treatment  based 
solely  on  the  molecular  weight  effect  would  predict  a  decreasing  effect  on  Cf  of 
mass  injection  with  increasing  molecular  weight  of  the  injected  gas.  This 
is  due  to  the  effect  of  the  molecular  weight  on  the  quantity  C.  However,  as 
can  be  seen  from  the  results  in  Tables  VI  and  VII,  the  present  investigation 
indicates  the  opposite  effect  with  regard  to  injection  of  argon,  and  carbon 
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dioxide.  This  result  is  due  to  the  fact  that  the  quantity  C  is  temperature 
dependent  and  the  heat  capacity  effect  of  the  injected  gas  must  be  taken  into 
account.  The  present  investigation  indicates  that  carbon  dioxide  injection 
has  a  greater  effect  on  reduction  of  c^  than  argon  injection.  Although  the 
molecular  weight  of  carbon  dioxide  is  greater,  the  temperature  distribution 
produced  by  its  injection  is  considerably  lower  than  that  for  argon  injection 
(see  Figures  8(b)  and  11(b))  due  to  its  high  heat  capacity;  hence,  the  net 
effect  on  Cf  is  greater  than  that  for  argon  injection.  It  is  expected,  however, 
that  the  heat  capacity  effect  would  decrease  with  decreasing  surface  tempera¬ 
tures  and  at  low  Mach  numbers  the  effect  of  injection  would  be  directly 
related  to  molecular  weight. 

In  addition  to  the  effect  of  including  TVC  on  the  skin  friction  predictions, 
interesting  effects  are  produced  regarding  predictions  of  heat  transfer.  The 
results  in  Table  VI  based  on  calculations  without  TVC  indicate  that  the  Stanton 
number  for  helium  injection  is  lower  than  that  for  any  of  the  other  gases.  For 
a  binary  mixture,  the  surface  heat  transfer  rate  (see  equation  (95))  is  com¬ 
posed  of  a  temperature  gradient  term  and  concentration  gradient  term.  Since 
the  heat  capacity  of  helium  is  much  greater  than  that  of  air,  the  concentration 
gradient  will  contribute  to  a  heat  flux  away  from  the  wall  and  the  temperature 
gradient  term  which  contributes  to  a  heat  flux  towards  the  surface  will  be 
lowered  by  the  heat  capacity  effect.  The  net  result  indicates  heat  transfer 
away  from  the  wall  for  the  set  of  conditions  investigated  here.  However, 
with  TVC  (see  Table  VII),  the  conduction  term  is  increased  to  the  point  where 
the  net  result  is  heat  transfer  to  the  wall.  Thus,  without  TVC,  helium  appears 
to  be  the  best  coolant  of  the  gases  studied  here,  and  with  TVC,  carbon  dioxide 
appears  to  be  the  best  coolant.  Again,  this  is  due  to  the  fact  that  the  effects  of 
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including  transverse  curvature  are  felt  the  most  with  helium  injection  and 
the  least  with  carbon  dioxide  injection. 
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TABLE  V 

COMPARISON  OF  COMPUTED  BOUNDARY  LAYER  PARAMETERS 
FOR  AIR  INJECTION  NEGLECTING  TYC  -  Pr  =  0.71 

P  v 

w  w  -  0.01434 


v  iS 

St  /cf 

**ll 

Present 

Method 

Reference  31 

Present 

Method 

Reference  31 

0.  20 

0.  5071 

0.  5088 

0.  0570 

0.  0581 

0.  45 

0.  5383 

0.  5327 

0.  0956 

0. 0962 

0. 60 

0.  5410 

0. 5490 

0.  1200 

0.  1308 

0.  72 

0. 5509 

0.  5573 

0. 1397 

0. 1508 

0.  86 

0. 5627 

0.  5675 

0. 1628 

0.  1737 

0.  90 

0. 5758 

0.  5705 

0. 1694 

0.  1801 

0.  96 

0. 5708 

0. 5752 

0.  1785 

0.  1897 

1.00 

0.  5740 

0.  5783 

0.  1846 

0. 1962 

TABLE  VI 

COMPUTED  BOUNDARY  LAYER  PARAMETERS, 
NEGLECTING  TVC  -  VARIABLE  FLUID  PROPERTIES 


He  0.3183  0.  0036  -0.1536  0.3695  0.0022  -0.2355 

Air  0.1398  0.0140  0.0747  0.1595  0.0114  0.0616 
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TABLE  VII 

COMPUTED  BOUNDARY  LAYER  PARAMETERS, 
INCLUDING  TVC  -  VARIABLE  FLUID  PROPERTIES 


P  v 
w  w 

~p~  u 
e  e 


0.  01434 


INJECTED 

GAS 

x/i  =  C.  75 

x/i  =  0.90 

6  *  / / 

cf 

St  x  10 

6  *  it 

Cf 

St  x  10 

He 

0.  1937 

0. 0147 

0.  6603 

0.  2475 

0. 0099 

0.  1656 

Air 

0. 1266 

0.  0232 

0.  1205 

0.  1434 

0.  0191 

0.  1001 

Ar 

0.  1280 

0. 0298 

0.  1723 

0.  1451 

0.  0249 

0.  1439 

co2 

0.  1060 

0. 0245 

0.  0529 

0.  1185 

0.  0203 

0.  0491 
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8.0  CONCLUDING  STATEMENTS 

The  method  described  in  this  report  has  proven  capable  of  solving  the 
binary  diffusion  laminar  boundary  layer  equations.  It  is  applicable  to  axi- 
symmetric  flows  with  or  without  transverse  curvature  effects  included  and 
to  two-dimensional  flows.  The  mathematical  transformations  used  make  it 
possible  to  easily  include  or  neglect  the  transverse  curvature  terms  within 
the  framework  of  one  computer  program.  A  combination  of  the  Probstein- 
Elliott  and  Levy-Lees  transformations  introduces  a  common  term  t  to  the 
system  of  equations  being  solved.  If  this  term  is  retained,  the  higher  order 
transverse  curvature  effects  are  accounted  for  in  the  numerical  solution  to 
the  equations;  if  this  term  is  neglected,  the  equations  assume  the  standard 
Levy- Lees  form  and  the  second-order  transverse  curvature  effects  are  not 
accounted  for.  The  general  method  of  solution  is  identical  for  both  cases; 
however,  the  proper  inverse  transformations  used  to  go  back  to  the  physical 
plane  depend  on  whether  or  not  the  quantity  t  is  included  or  neglected  in  the 
solution  of  the  partial  differential  equations. 

The  method  of  solution  has  been  programmed  for  the  IBM  7094  com¬ 
puter  for  the  following  mixtures:  helium-air,  argon-air,  and  carbon 
dioxide-air.  The  program  also  contains  an  option  which  bypasses  the 
solution  of  the  species  equation  when  a  calculation  with  pure  air  is  desired. 
Variable  fluid  properties  given  as  functions  of  local  temperature  and  mass 
fraction  throughout  the  boundary  layer  are  used.  Although  the  equations 
appear  to  be  more  complex  when  second-order  transverse  curvature  is 
accounted  for,  the  actual  computing  times  are  less.  This  is  due  to  the  fact 
that  convergence  of  trial  solutions  of  the  momentum  equation  is  more  rapid 
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when  the  transverse  curvature  term  appears  in  the  equation.  Approximately 
21  stations  along  the  cone  were  used  to  obtain  solutions  neglecting  the  TVC 
term  and  27  stations  were  used  when  it  was  included.  The  total  computing 
times  for  the  complete  body  were  about  47  minutes  and  36  minutes  respectively. 

The  results  indicate  that  the  TVC  effect  is  significant  and  increases 
with  decreasing  molecular  weight.  As  a  corollary  of  this,  the  effect  of 
helium  injection  on  the  reduction  of  skin  friction  relative  to  the  effect  of  air 
injection  is  far  less  pronounced  when  TVC  is  taken  into  account.  Solutions 
to  the  equations  containing  the  TVC  term  yield  values  of  skin  friction 
coefficient  and  Stanton  number  higher  and  values  of  displacement  thickness 
lower  than  solutions  which  neglect  this  term. 

The  results  obtained  indicate  the  necessity  of  using  exact  fluid 
properties  to  ascertain  the  relative  effects  of  injection  of  various  gases. 

The  investigation  of  viscous  interaction  for  hypersonic  flow  over 
slender  cones  with  mass  injection  is  recommended  for  future  study.  This 
investigation  could  be  conducted  by  using  the  present  method  together  with 
an  iteration  between  external  pressure  distribution  and  computed  displace¬ 
ment  thickness.  It  is  felt  that  the  viscous  interaction  effect  will  increase 
with  decreasing  molecular  weight  of  the  injected  gas  and,  in  general,  will 
cause  behavior  similar  to  that  caused  by  inclusion  of  the  transverse 
curvature  effect. 
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Figure  5. 
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Figure  8.  -  Continued,  (b)  Temperature  Profiles 
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Figure  8.-  Concluded,  (c)  Mass  Fraction  Profiles 
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Figure  11.  -  Comparison  of  Boundary  Layer  Profiles  at  x//  =  0.  75, 
Including  TVC.  (a)  Dimensionless  Velocity  Profiles 


80 


AEDC-TR-46-183 


T 


2400 

2000 

1600 

(°R) 

1200 

800 

400 


0  20.0  40.0  60.0  80.0 


Y 


Figure  11.  -  Continued,  (b)  Temperature  Profiles 
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Figure  li.  -  Concluded,  (c)  Mass  Fraction  Profiles 
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